Dilepton to photon ratio, a viscometer of QGP 
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In the Israel-Stewart's 2nd order hydrodynamics, a viscous effect on dilepton emission from a 
QGP medium is investigated. Dileptons are strongly affected by QGP viscosity. Large invariant 
mass dileptons, due to their lower velocity, are less affected by viscosity than the low invariant mass 
dileptons. We also show that the ratio of photon to dilepton is sensitive to the viscosity and can 
serve as a viscometer for QGP. 
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I. INTRODUCTION 

It is expected that collisions between two nuclei at 
ultra-relativistic energies will lead to a phase transition 
from hadrons to the fundamental constituents, quarks 
and gluons, usually referred to as Quark-Gluon-Plasma 
(QGP). Although the nature of the phase transition still 
remains largely uncertain, it is an acceptable reality that 
QGP is formed. 

Electromagnetic signals, photons and dileptons are ex- 
cellent probe of QGP, formed immediately after the col- 
lision of the two nuclei. The interaction, being primarily 
electromagnetic, photons and dileptons tend to retain the 
information of early times rather more efficiently com- 
pared to hadronic signals; the hadrons interact strongly 
and thus tend to erase the information of early time re- 
taining the information only around the freeze-out sur- 
face. 

In ideal hydrodynamic models, photon and dilepton 
production in high energy nuclear energy has been stud- 
ied extensively However, it is now realized that 
the strongly interacting medium, produced in Au+Au 
collisions, must be treated as a viscous medium. Gravity 
dual (AdS/CFT) theories suggest that specific viscosity, 
i.e. viscosity to entropy ratio of any matter has a lower 
bound, the so called KSS bound 77/s = l/4n [alH- Even 
though, photons and dileptons are important probe of 
QGP matter, viscous effects on these signals still need to 
be investigated extensively. Dusling [3, |8| studied viscous 
effects on these signals. In a recent paper [9], we have in- 
vestigated viscous effects on photons produced from the 
QGP medium. It was shown that viscous effects on pho- 
tons are strong. It was suggested that the slope parame- 
ter of the transverse momentum spectra for photons can 
possibly limit the thermalisation time and viscosity. In 
the present paper, we investigate viscous effect on dilep- 
ton emission from QGP. Unlike photons, dileptons are 
quite massive. Invariant mass of the dileptons make them 
unique probe of QGP. In addition to transverse momen- 



tum, the invariant mass of the dileptons can be tuned to 
explore very early QGP. For example, dileptons of invari- 
ant mass Mi nv ~ 1 GeV are expected to be formed in the 
time scale T/ or ~ l/Mi nv ~ 0.2 fm and can explore the 
QGP at very early time r ~ 0.2 fm. Viscous effects on 
dileptons will also depend on the invariant mass. Viscous 
drag on dileptons is proportional to relative velocity be- 
tween the fluid and the lepton pairs. Depending on the 
invariant mass, viscous drag will vary. As will be shown 
later, this unique feature makes them very sensitive to 
the viscosity. Interestingly enough, we find in this work 
that the ratio, dileptons to photons is quite sensitive to 
the viscosity and thus the ratio turns out to be an excel- 
lent viscometer. 



II. HYDRODYNAMICAL EQUATIONS, 
EQUATION OF STATE AND INITIAL 
CONDITIONS 

In a hydrodynamical model, the invariant distribution 
of dileptons is obtained by convoluting the dileptons pro- 
duction rate with space-time evolution of the fluid. We 
assume that in ^/s NN —200 GeV, Au+Au collisions at 
RHIC, a baryon free QGP fluid is formed. Space-time 
evolution of the fluid is obtained by solving 2nd order 
Israel-Stewart's theory, 
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Eq[T] is the conservation equation for the energy- 
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momentum tensor, T^ v = (e + p)u^u v — pg 1 
e, p and u being the energy density, pressure and fluid 
velocity respectively. is the shear stress tensor 

(we are neglecting bulk viscosity). Eqj2] is the relax- 
ation equation for the shear stress tensor . In Eq|3J 
D = u^dfj, is the convective time derivative, V <M u ly> — 
^(V^u y + V u u^) - \{d.u){g^ - u"u") is a symmetric 
traceless tensor, r\ is the shear viscosity and tv is the 
relaxation time. It may be mentioned that in a confor- 
mally symmetric fluid relaxation equation can contain 
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additional terms ■ Assuming boost- invariance, EqsQ] 
and [5] are solved in (r = \Jt 2 — z 2 , x, y, t] s — jlnf^f) 
coordinates, with the code " 'AZHYDRO-KOLKATA" ', 
developed at the Cyclotron Centre, Kolkata. Details of 
the code can be found in 

Eqs ll|2l are closed with an equation of state p = 
p(e). Lattice simulations |L4T - [l7| indicate that the 
confinement-deconfinement transition is a cross over, 
rather than a 1st or 2nd order phase transition. There is 
no critical temperature for a cross-over transition. The 
inflection point on the Polyakov loop is generally de- 
fined as the pseudo critical temperature. In Wuppertal- 
Budapest simulation [l|| [13], the pseudo critical tem- 
perature T c «170 MeV. In the present simulation, for 
the QGP phase, we use an EOS based on Wuppertal- 
Budapest simulation. 

Solution of partial differential equations fEqs lll2[) re- 
quires initial conditions, e.g. transverse profile of the 
energy density (e(x,y)), fluid velocity (v x (x, y), v y (x, y)) 
and shear stress tensor (tt^ (x,y)) at the initial time Tj. 
One also need to specify the viscosity (77) and the relax- 
ation time (r w ). A freeze-out temperature is also needed. 
In the following, we will study viscous effects on dilepton 
production from the QGP phase only. The hydrodynam- 
ical equations are then solved till the freeze-out temper- 
ature Tf = T c =170 MeV. At the initial time n, initial 
energy density is assumed to be distributed as [18| 



where b is the impact parameter of the collision. N part 
and N co u are the transverse profile of the average partic- 
ipant and collision number respectively, calculated in a 
Glauber model. fhard=0.13 is the hard scattering frac- 
tion [19]. £i is the central energy density of the fluid in 
impact parameter 6 = collision. In the following, we as- 
sume that the fluid is thermalised at the time scale Ti=0.6 
fm to central temperature To=350 MeV. A large number 
of charged particle's data e.g. identified particles spec- 
tra, elliptic flow etc. data are explained in hydrodynami- 
cal model with similar initial time and temperature scale 
fl8j |. We also assume that initial fluid velocity is zero, 
v x (x,y) = v y {x,y) — 0. The shear stress tensor was ini- 
tialized with boost- invariant value, tt xx — ir vv — 2r]/3Ti, 
tt xv =0. For the relaxation time, we use the Boltzmann 
estimate tv = Sri/Ap. We further assume that shear vis- 
cosity to entropy ratio is a constant throughout the evo- 
lution and we simulate Au+Au collisions for a range of 
■q/s, r?/s=0-0.12. 



III. DILEPTON PRODUCTION RATES 



e(h,x, y) = Ei[{l-f har d)N part (b,x,y)+f hard N co u(b,x,y)], 

(3) 



The rate of dilepton emission, from a QGP, due to 
quark-antiquark annihilation can be written as, 
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where q = (qo,q) is the virtual photon 4- momenta and function can be approximated as, 
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(pi 2 + P2 ) is the invariant mass, 
is the relative velocity between 



In EqU v 12 = 

the quark-antiquark pair and a = — -g^f? — - is the qq 
cross-section. f(E,T) in EqU is the quark/antiquark 
distribution function. In ideal fluid evolution, QGP is in 
(local) equilibrium and f(E,T) = f eq (E,T) = 1+ ^ E / T ■ 
For E/T — p.u/T << 1 the distribution function can be 
approximated by a Boltzmann distribution and integra- 
tions in EqU can be analytically performed, 
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In 



f eq {l + 5f neq ), Sf neq « 1. 



(6) 



The non-equilibrium correction 5f neq depend on dissi- 
pative forces as well as on particle momenta. For shear 
viscosity, non-equilibrium correction to equilibrium dis- 
tribution can be obtained from the following ansatz, 
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Dilepton production rate from non-equilibrium QGP 
then has two parts, an equilibrium part and a non- 
equilibrium part, 



Unlike an ideal fluid, in viscous evolution, distribution 
functions / is modified due to non-equilibrium correc- 
tion. For small non-equilibrium effect, the distribution 
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FIG. 1: (color online) Transverse momentum spectra for 
dileptons of invariant mass M=600 MeV, from evolution of 
ideal and minimally viscous (?7/s=0.08) QGP. For viscous 
fluid, spectra obtained with and without the non-equilibrium 
correction to distribution function, are shown separately. In 
the bottom panel, ratio of production from viscous fluid with 
and without the non-equilibrium correction to the ideal fluid 
is shown. 



Approximating the equilibrium distribution function 
by Boltzmann distribution, Dusling and Lin [7( evaluated 
the non-equilibrium correction to the dilepton production 
rate, 

\dSJ neq 3(e + p)T2 W 

Since non-equilibrium correction to distribution func- 
tion is assumed to be small, it is essential that, 

(wL <K (l0 > 



IV. EFFECT OF VISCOSITY ON DILEPTON 
SPECTRA 

For viscosity to entropy ratio r//s=0, 0.04, 0.08 and 
0.12, we have simulated 20-40% Au+Au collisions and 
computed dileptons invariant distribution. In FigJTJ 
transverse momentum spectra for dileptons of invariant 
mass M=600 MeV, from evolution of ideal and mini- 
mally viscous (?7/s=0.08) QGP fluid, are shown. In vis- 
cous evolution, dilepton spectra is modified due to (i) 
changed space-time evolution and (ii) non-equilibrium 
correction to the equilibrium distribution function. To 



FIG. 2: (color online) Transverse momentum spectra for 
dileptons of invariant mass M=300 MeV, form evolution of 
QGP fluid with viscosity to entropy ratio 7//s=0, 0.04, 0.08 
and 0.12. In the lower panel, ratio of non-equilibrium correc- 
tion to equilibrium production is shown 

understand the effect of viscosity, spectra, with or with- 
out the non-equilibrium correction to the equilibrium dis- 
tribution function are shown separately. In the lower 
panel of FigJTJ the ratio dileptons from viscous QGP and 
ideal QGP are shown. If the non-equilibrium correction 
to the equilibrium distribution function is neglected, in 
viscous evolution dilepton yield increase by a factor of 
^1.3 . The increase is largely pt independent. The ratio 
become px dependent when the non-equilibrium correc- 
tion is included. It is expected that the non-equilibrium 
correction increases quadratically with p? (see EqQ. 

Simulated transverse momentum spectra for dileptons 
of invariant mass M=300 MeV, for four values of vis- 
cosity to entropy ratio, rj/s— (ideal fluid), 0.04, 0.08 
(AdS/CFT lower bound) and 0.12 are shown in FigH 
In viscous evolution, dilepton spectra is stiffened, stiff- 
ening is more for more viscous fluid. The result is un- 
derstood. In viscous evolution, non-equilibrium correc- 
tion introduces additional momentum dependence which 
is linearly proportional to shear stress tensors. One also 
note that at very low px, yield is marginally reduced in 
viscous fluid. This result is not unphysical. At low px 
non-equilibrium correction to the equilibrium distribu- 
tion function can be negative. A similar decrease is seen 
in case of low pt hadron production also [2(| ■ 

As it was mentioned earlier, it is essential that the non- 
equilibrium correction to dilepton production is small 
compared to the equilibrium production. The non- 
equilibrium distribution function is defined under the as- 
sumption that the system is not far away from equilib- 
rium and the non-equilibrium correction to the equilib- 
rium distribution function is small {f neq << 1). In the 
lower panel of FigJ5J the ratio of non-equilibrium correc- 
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FIG. 3: (color online) same as in Fig[2] but for dilcpton in- 
variant mass M=600 MeV. 



tion to dilepton production to the equilibrium production 
is shown. At very low px , px < OSMeV, non-equilibrium 
correction is negative but small. At larger px correction 
is positive and rapidly increase with px- For viscosity to 
entropy ratio 77/3=0. 04, 0.08 and 0.12, non-equilibrium 
correction equals the equilibrium production at px ~ 2.9, 
1.9 and 1.5 GeV respectively. Definitely, viscous hydro- 
dynamics is inapplicable beyond those px range. Indeed, 
if we assume that hydrodynamics remain applicable till 
\SN nec t/N £ i\ < 0.5, the dilepton production from QGP, 
with viscosity to entropy ratio, 0.04, 0.08 and 0.12 can 
be computed with confidence only up top-r=(0-2.1 GeV), 
(0-1.5 GeV) and (0-1.2 GeV). The result is similar to that 
obtained in our simulations for photon production in vis- 
cous hydrodynamics @. Viscous effects are also very 
strong on photon and photon production can be com- 
puted only in a limited px range. 

In Fig|3] same results are shown for dilepton invariant 
mass M=600 MeV. As before, at low px, non-equilibrium 
correction contribute negatively, and dilepton production 
is reduced in viscous evolution. The correction is positive 
at higher px and dilepton yield is increased. In the lower 
panel, the ratio of non-equilibrium correction to equilib- 
rium production is shown. As expected, correction in- 
creases with px and also with viscosity. For viscosity to 
entropy ratio rj/s—Q. 04, 0.08 and 0.12, non-equilibrium 
corrections become equal the equilibrium contribution at 
px ~ 3 GeV, 2.1 GeV and 1.7 GeV. We have not shown 
but similar result is obtained for higher mass M=900 
MeV dileptons. For M=900 MeV, for viscosity to entropy 
ratio ?7/s=0.04, 0.08 and 0.12, non-equilibrium correction 
equals the equilibrium contribution at px ~ 3.2 GeV, 
2.4 GeV and 2 GeV. It is interesting to note that the 
transverse momentum px ma!C at which non-equilibrium 
correction equals to equilibrium correction, in addition 
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FIG. 4: (color online) elliptic flow in ideal and viscous QGP. 
For viscous QGP, flow with and without the non-equilibrium 
correction is shown separately. 



to viscosity, also depend on the dilepton invariant mass. 
For example, for ?y/s=0.12, non-equilibrium contribution 
equals the equilibrium contribution at Px mal — 1.5 GeV, 
1.7 GeV and 2 GeV respectively for dilepton invariant 
mass M=300, 600 and 900 MeV. It appear that at large 
px viscous effect is lessened in more massive dileptons. 
The result can be understood qualitatively. Dilepton ve- 
locity is approximated as 1/ y/ M 2 + p T . More massive 
dileptons will move slowly, consequently, viscous drag on 
massive dileptons is less than on lighter dileptons. The 
feature is important. It will be shown later this partic- 
ular feature can be utilised to accurately measure QGP 
viscosity. 



V. EFFECT OF VISCOSITY ON DILEPTONS 
ELLIPTIC FLOW 

Elliptic flow is an important observable in relativistic 
energy nuclear collisions. Large elliptic flow of observed 
hadrons is indirect proof of thermalisation of the system. 
In FigHl we have shown the simulation results for el- 
liptic flow for dileptons of invariant mass M=600 MeV. 
Results are shown for ideal and viscous ?7/,s=0.08 evo- 
lution. Flow with and without the non-equilibrium cor- 
rection is shown separately to understand the effect of 
viscosity. In ideal fluid evolution, elliptic flow increases 
with increasing px up to px ~2.5 GeV. At larger px, 
flow decreases, indicating that highp^ dileptons are over- 
whelmingly from the early non-flow phase. Elliptic flow 
dileptons is not large, even in ideal QGP, in 20-40% col- 
lision, the highest V2(px — 2.5GeV) rj 0.04. Interest- 
ingly, if the non-equilibrium correction is neglected, in 
viscous fluid evolution, elliptic flow increases. However, if 
the non-equilibrium correction is included, the simulation 
correctly predicts reduced flow. The result indicates the 
importance to have a consistent theory. Entirely wrong 
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FIG. 5: (color online) Elliptic flow for dileptons of invariant 
mass M=600 MeV, as a function of QGP viscosity. 
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FIG. 6: For dileptons of invariant mass M=900 MeV, the 
ratio of dileptons to photon is shown. The four lines (top 
to bottom) corresponds to viscosity to entropy ratio rj/s=0, 
0.04, 0.08 and 0.12 respectively. 



conclusion can be reached in an inconsistent theory. 

In FigJSl elliptic flow dileptons of invariant mass 
M=600 MeV is shown as a function of viscosity. The 
flow is shown only in the pr range when viscous hydro- 
dynamics is applicable. With increasing viscosity flow is 
reduced and become negative at large viscosity to entropy 
ratio T]/s=0.12. 



VI. DILEPTON TO PHOTON RATIO AS 
VISCOMETER FOR QGP 

QGP viscosity is an important parameter. Theoretical 
estimates of the ratio, (shear) viscosity over the entropy 
density, i]/s cover a wide range, 0.0-1.0. String theory 
based models (ADS/CFT) give a lower bound on vis- 
cosity of any matter n Is > 1/4tt |5j. In a perturbative 
QCD, Arnold et al ^estimated rj/s ~ 1. In a SU(3) 
gauge theory, Meyer |22| gave the upper bound r]/s <1.0, 
and his best estimate is ?7/s=0.134(33]at T = 1.165T C - 
At RHIC region, Nakamura and Sakai [23| estimated the 
viscosity of a hot gluon gas as 7//s=0. 1-0.4. Attempts 
have been made to estimate QGP viscosity directly from 
experimental data. Gavin and Abdel-Aziz [24| proposed 
to measure viscosity from transverse momentum fluctua- 
tions. From the existing data on Au+Au collisions, they 
estimated QGP viscosity as rj/s—Q. 08-0. 30. Experimen- 
tal data on elliptic flow has also been used to estimate 
QGP viscosity. Elliptic flow scales with eccentricity. De- 
parture from the scaling can be understood as due to 
off-equilibrium effect and utilised to estimate viscosity , 
fy/s=0.11-0.19 [H-ii]. The estimates are well within the 
upper bound obtained rj/s < 0.5 obtain in [29ll30l|. How- 
ever, hydrodynamical models estimates of QGP viscosity 
can be largely uncertain due to uncertainty in initial con- 
ditions. In a detailed analysis of RHIC data on phi meson 
11,], it was shown that uncertainties initial conditions in 
hydrodynamic models can induce very large uncertainty 



(~175% or more) in viscosity to entropy ratio. 

Ideally speaking, hadrons are not the best probe to 
study QGP viscosity. Hadrons arc emitted from the 
hadronic phase. They can carry viscosity information of 
the QGP phase only indirectly that too with the assump- 
tion that the non-perturbative hadronisation process did 
not erase the memory of the earlier QGP phase. Elec- 
tromagnetic probes like photons and dileptons, on the 
other hand are emitted directly from the QGP phase 
and can carry directly the information of viscosity of 
the QGP phase. However, below the critical tempera- 
ture, QGP fluid transforms into hadronic fluid. Pho- 
tons and dileptons are also produced from the hadronic 
phase. Experimental data then contain contributions of 
both the QGP and hadronic phases. Only in a limited 
phase space QGP phase dominates the production. For 
photons, QGP phase dominates the production in the pr 
range, 1.5 < pr < 2.5 GeV 0, BJ. Large invariant mass 
dileptons are also dominated by the QGP phase. 

In a recent paper Q , we have studied viscous effect on 
direct photon production. Details of the calculation can 
be seen in [9[ . We now argue that the dilepton to photon 
ratio can serve as an effective viscometer for QGP. The 
argument is simple, spectra of both photon and dilep- 
tons are stiffened in viscous QGP. Stiffening is more in 
more viscous QGP. Low invariant mass dileptons behaves 
essentially as photons and are similarly affected by the 
viscosity. Viscous effects will largely cancels out and the 
ratio will not be sensitive to viscosity. But large invari- 
ant mass dileptons moves slowly, viscous drag on them 
is less than on photons. Photon spectra will stiffen more 
than the large invariant mass dilepton spectra and the 
ratio will drop with increasing viscosity. The ratio has 
the advantage that the uncertainties in initial conditions 
can be largely eliminated. 

In FigJSl for viscosity to entropy ratio r//s=0, 0.04, 0.08 
and 0.12, the ratio of dilepton to photon production, for 
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FIG. 7: For dileptons invariant mass M=300, 600 and 900 
MeV, dileptons to photon ratio, at pr=l-5 GeV, is shown as 
a function of viscosity. 

dileptons invariant mass M — 900 MeV, is shown. Wc 
have shown the results only in the pr range in which 
viscous hydrodynamics remain applicable for the high- 
est r]/s=0.12 considered here. The ratio depends weakly 
on pt, pt dependence of the non-equilibrium correction 
largely canceling out. As argued before, the ratio drops 
with increasing viscosity. The dependence of viscosity 
can be seen more clearly in Fig{7] In FigEl the dilepton 
to photon ratio at pr ~1.5 GeV is shown as a function 
of viscosity. For invariant mass M=300 MeV, the ra- 
tio hardly changes with of viscosity. For invariant mass 
M=600 MeV, the ratio decreases with viscosity, the de- 
crease is even faster for invariant mass M=900 MeV. Fast 
change in the ratio for large invariant mass dileptons 
makes it very sensitive to viscosity. For dilepton mass 
M=600 MeV, if the ratio is measured within 10% accu- 
racy, viscosity to entropy ratio can be estimated within 
an accuracy of ~5%. The sensitivity is increased to ~2% 



for invariant mass M=900 MeV. We have shown the ratio 
at a fixed pt ~ 1-5 GeV. The slope of the ratio is largely 
Pt independent and sensitivity of the ratio to viscosity 
remain unaltered at other pt also. It may be mentioned 
here that the above discussion is based on the assumption 
that pt ~1.5 GeV, photons and dielptons are dominantly 
from the QGP phase. While ideal hydrodynamics anal- 
ysis support the assumption, it quite likely that window 
for QGP may well get modified with viscous effect. 

It has been pointed out before [3l|, [32| that the dilep- 
ton to photon ratio is a clear signal of QGP; by taking 
the ratio uncertainties of initial conditions and other pa- 
rameters largely get cancelled out. Now, with the ratio 
being an accurate viscometer, the advantage of the ratio 
can be further utilised. 

VII. SUMMARY 



To summarize, we have studied viscous effects on dilep- 
ton production from QGP. In viscous dynamics, dilepton 
production is modified due to (i) changed space-time evo- 
lution of the fluid and (ii) non-equilibrium correction to 
the equilibrium distribution function. With initial condi- 
tions appropriate for -^5^^=200 GeV Au+Au collisions, 
space-time evolution of QGP was obtained by solving 
Israel-Stewart's 2nd order hydrodynamics. Invariant dis- 
tribution for dileptons was obtained by convoluting the 
dilepton production rates over the space-time evolution. 
Effect of viscosity to stiffen the dilepton spectra and re- 
duce elliptic flow. Stiffening is more in more viscous 
QGP. Elliptic flow is also reduced more in more viscous 
fluid. It was also shown that viscous effect on dileptons 
is strong. Even for minimally viscous QGP (j]/s — 0.08), 
dilepton production can be reliably computed, only in a 
limited pt range. It is also indicated that for large mass 
dileptons, the ratio of dilepton to photon yield is sensi- 
tive to viscosity to entropy ratio. Accurate estimate of 
QGP viscosity can be obtained by measuring the ratio 
experimentally. 
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